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SUMMARY 


The  error  function,  which  occurs  in  much  of  the  literature  of  mathemeticn,  phyaica,  and 
engineering,  ia  generalised  to  handle  the  normal  probability  distribution  in  n  dimensions. 

Explicit  integral  repreaentations  for  theae  functions  are  found  to  be  of  two  general  forma, 
depending  upon  whether  n  ia  even  or  odd. 

Some  readily  eatablished  recursion  formulas  and  other  relationahipa  are  derived  for  theae 
functions. 
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The  mathematical  theory  of  probability  and  the  related  techniques  of  statistics  are  being  used 
by  an  increasing  number  of  workers  in  many  diverse  fields  embracing  the  sciences  and  engineering, 
as  well  as  mathematics. 

Of  particular  importance,  therefore,  to  engineers,  to  operations  and  systems  analysts,  and  to  the 
designers  of  experiments  are  certain  standard  probability  distributions.  The  most  widely  employed 
of  the  continuous  distributions  is  undoubtedly  the  “Gaussian,"  or  “normal"  distribution,  which  is 
of  enormous  theoretical,  historical,  and  practical  importance. 

The  normal  distribution,  with  zero  mean,  is  given  by 

p  (x)  = _ L_  e~~2°1  (1) 

/27a 

and  is  completely  specified  by  the  parameter  a,  called  the  “standard  deviation." 

This  distribution  has  been  generalized  to  n  dimensions.  The  most  general  n -dimensional  nor¬ 
mal  distribution  contains  parameters  to  account  for  nonzero  means  for  the  n  independent  variables, 
for  correlations  among  the  variables,  and  for  unequal  standard  deviations  with  respect  to  each  of 
the  n  variables. 

This  report  will  concern  itself  with  a  very  special  case  of  the  n-dimensional  distribution.  In 
particular,  the  means  will  all  be  assumed  zero,  the  correlations  will  all  be  assumed  zero,  and  the 
standard  deviations  will  be  assumed  equal.  The  resulting  probability  distribution  is  then  given  by 


1  -  *1  •••  *n 

p(Xj.  *2 . xn)  -  - ^ - e~  2„i  (2) 

(2  n)7«n 

In  the  one-dimensional  case,  a  commonly  occurring  expression  containing  p  (x)  is  that  for  the 
probability  with  which 

I  x  |  <  a 

This  expression  is  (for  a_-  0) 


(3) 

(3) 

(since  the  integrand  is  even) 


3 


This  expression  cannot  be  evaluated  in  closed  form  for  an  arbitrary  upper  limit,  but  occurs  so 
frequently  that  its  values  have  been  tabulated,  and  the  expression  itself  has  been  given  the  name 
of  “error  function.” 


The  customary  definition  of  the  error  function  is 


erf  (x)  s  — §.  f  e'y1  dy 

In  terms  of  this  definition,  it  follows  that 
Prob  1  |  x  |  <  A  I  ■= 


(4) 


(5) 


The  analog  of  the  problem  in  n  dimensions  presents  no  new  problems  in  rectangular  coordinates, 


since 


Prob 


*1  I  i  i  l  x2  I  — 


as  a  consequence  of  the  appropriate  integral  separating  into  a  product  of  integrals  with  respect  to 
one  variable. 


Something  new  does  arise,  however,  from  regarding  r  =  (xj,  ...»  xn)  as  a  vector  in  n  dimen¬ 
sions,  and  asking  what  is 

Prob  I  |  r  ]  v  A  i 

This  is  a  very  natural  question  to  ask  for  the  cases  n  -  2,  3. 

It  is  the  purpose  of  the  present  paper  to  investigate  this  question.  The  result,  as  will  be  seen, 
is  to  define  an  “error  function”  generalized  to  n  dimensions,  in  terms  of  which  the  required  proba¬ 
bility  can  be  written  in  a  manner  analogous  to  equation  (5). 

Some  simple  properties  of  these  generalized  error  functions  are  proved,  and  graphs  of  these 
functions  are  presented. 


« 


THE  SPECIAL  CASE  OF  TWO  DIMENSIONS 

As  a  natural  way  of  introducing  generalized  error  functions  it  is  instructive  to  consider  the 
case  of  two  dimensions.  This  case  is  well  known  among  the  users  of  probability  theory  because 
of  its  frequent  occurrence  and  the  fact  that  the  mathematics  fortuitously  permits  a  solution  in 
closed  form. 

The  two-dimensional  treatment,  moreover,  is  capable  of  a  direct  generalization  to  n  dimen¬ 
sions,  as  will  be  seen  later,  and  it  is  therefore  profitable  to  dwell  at  some  length  upon  this 
special  case.  This  will  now  be  done. 
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Introduce  polar  coordinates  (r,  &).  Then 


Prob  ||  r  |  <  a  |  -//i  e  2cr2  dx  dy 


2  no* 

clrcl*  about 
origin,  rodlut  & 
-2  7T 


-/V'jL.* 

T)  •'o  2rra2 


C  2a* 


r  dr  d  0 


( 


—  e  j  cr  2  r  dr 

a2 


(7) 


.•.Prob  I  |  r  |  <  a  I  =  1  -e“3<72 


(8) 


This  is  a  well  known  expression  which  appears  very  frequently  in  the  literature. 
The  above  result  suggests  defining  an  error  function  in  two  dimensions  by 

erf2  (x)  r  1  -e"* 

Then 
Prob  I  | 


(9) 


(10) 


in  two  dimensions,  analogously  to  equation  (5). 

The  above  computational  procedure  will  now  be  extended  to  an  n-dimensional  distribution.  It 
will  be  seen  later  that  all  the  even-dimensional  error  functions  are  expressible  in  closed  i'orm, 
although  not  always  conveniently  so. 


THE  GENERAL  CASE  OF  n  DIMENSIONS 

The  case  of  n  dimensions,  for  n  >  2,  is  carried  out  analogously  to  the  above  two-dimensional 
treatment. 


It  is  necessary  to  evaluate  the  following  expression: 


Prob  \  \  t  \  <  Ui  -J J  f - L - 

hyp»r«phere  (2  n  )  2  „  n 
■bout  origin. 


dx. 


dxn 


(ID 


The  integration  is  simplified  by  introducing  hyperspherical  coordinates,  (R,  0 j,  02» 
according  to  the  transformation  equations 


* 


-  R  ain  am  02 

X2  “  R  am  flirt  ^2 

|  =  R  ain  flj  ain  Oj 


n-2  n-1 


n-2 


n— 1 


a  in  co  a 


n-2 


-  R  a  in  tlj  a  in  02  •  •  •  ain  °»-4  cod  ^n-3 


(12) 


Xn->2  —  R  am  flj  caa 


X  =  R  COfl  ^ 


Equation  (11)  then  becomes 


Prob  t  I  r  I  <  a  i  •- 


1 


(2 n)1  an  0 


f  e'^Rn'ldR/  /  d[1-„ 


where  d  lln  is  an  element  of  hypersolid  angle  and  is  independent  of  R  and  d  R. 
In  particular. 


(13) 


CD  j 

1  =  Prob  I  |  r  |  <  u  U - L -  f  e’^Rn’‘dR f  f  d  fln 

(2n)^a»  0  ... 


(14) 


It  follows  that 


/  J‘VjT^ 

/  e—R"-1 

•'a 


dR 


Substituting  into  equation  (13), 


J  ~  if 

Prob  1  |  r  |  <  a  l  =  - j 


C4  2 

e*rb  Rn"‘  dR 


dR 


An  obvious  change  of  variable  in  the  integrals  gives  the  result 


Prob  I  l  r’  |  •_  1 1  -„r- 


r’yp  2  , 

J  e-u  u"'1  du 


L 


e”u  u-'du 


i 


(15) 


(16) 


(17) 

I 
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Now  define  a  generalized  n-dimensionai  error  function,  erfn(x),  by 


erf(x) 


f  e"“J  u11"1  du 

■4 

:  — « - - - 

f  e”u  u”"1  du 


(18) 


Equations  (17)  and  (18)  then  imply  that 
Prob  I  |  r  |  <  a  1  =  erfn  / a  \ 

WV 

which  is  analogous  to  equation  (S),  and  is  valid  for  n  =  1,  2,  3,  4,  .  .  . 

For  n  =  1,  2,  the  definition  (18)  reduces  to  the  definitions  (4)  and  (9),  respectively. 

SOME  FURTHER  ANALYSIS 

Consider  the  integral 


(19) 


CC 

Ik  3  J  e'“J  uk  du 


(20) 


which  is  used  as  a  normalizing  factor  in  equation  (18). 
Integrating  by  parts,  it  is  readily  established  that 


i  -JllJL  I 
*k  *2 - 2  *k 


It  follows  by  mathematical  induction  that 
Ijju  ~  m  !  Ij  #  m  -  0,  1,  2,  ... 


(21) 


(22) 


and 


;  m  - 0.  1.  2 - 

22m  m! 


(23) 


but 


I,=  f  e”u  u  d  u  =  -L 

1  L  2 


(24) 


and 


l0  =/  e““2du 


(25) 
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so  that 


=3m! 


and 


t  (2  m)!  Jn 
2nl  ~  22m+1  m  ! 

The  definition  (18)  can  therefore  re  be  written 


(26) 


(27) 


erf2mnW=- 


)2m 


y^”(2m)! 


u2m  d u  ;  m  =  0,  1,  2,  .  .  . 


(28) 


and 


erf2m(X)  = 


! _ [*  e"“J  u2  S2ni“1  du  ; 

1  )  !  Jn 


(m  -1  ) !  Jq 


m  ~  1^  2,  3,  .  .  ■ 


(29) 


Equations  (28)  and  (29)  may,  ift  A  desired,  together  be  taken  as  the  definition  of  the  generalized 
error  function,  rather  than  equation  x  (18). 

The  error  functions  will  now  ba«»e  shown  to  satisfy  certain  relationships. 

Integrate  equation  (28)  by  parts^ts,  The  resulting  expression  simplifies  to 


erf2mnW  = 


(2x)2m-m!e^;+erf2m+3(x)  ;  .  ,  0,  1,  2,  .  .  . 


yiT(2m  +  1)! 

Proceeding  similarly  with  equalisation  (29),  one  gets 


(30) 


erf2m<X>  - 


x2me"' 


f  erf. 


*2m»2£2^  ;  m  -  1,  2,  3,  .  .  . 

Recalling  equation  (9),  it  folio  .ows  by  mathematical  induction  on  equation  (31)  that 


(31) 


erf2m(x)  -  1  -  e' 


,2 

Ti 


2! 


1  +2LI  +  .  .4  x 


2(m-l) 


(m-1) 


;  m  0,  1,  2,  .  .  . 


(32) 


Equation  (32)  shows  that  all  t_»  the  even-dimensional  error  functions  are  expressible  in  closed 
form,  although  for  sufficiently  largjjge  m  the  closed  form  expressions  become  increasingly  compli¬ 
cated. 


It  follows  also  from  equation  •  (32)  that 

2 

1  -e"  erfjfx)  j  erf4(x)  -  •  erf6(x)  >  .  .  . 

and  that  for  any  preassigned  x , 


(33) 
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In  s  similar  manner,  mathematical  induction  applied  to  equation  (30)  gives  the  result 


erf2m+l(X>  =  etfto-^r 


(2x)  0  !  j.  (2x)3  1 !  ^  (2x)2m-5(m-l)!  . 

l!  3!  *’*  (2m-l) ! 


(35) 


m  =  1,  2,  3,  .  .  . 


Since  erft(x)  =  erf(x)  cannot  be  expressed  in  closed  form,  neither  can  the  odd-dimensional 
error  functions. 


However,  it  is  clear  from  equation  (35)  that 

erf(x)  *  erfj(x)  >  erf3(x)  >  erf5(x)  >  .  .  .  (36) 


since  x  was  assumed  non-negative  in  the  definition,  equation  (18),  in  view  of  the  probability 
application  with  regard  to  which  the  generalized  error  functions  were  introduced. 

From  the  point  of  view  of  the  pure  mathematician,  it  is,  of  course,  desirable  to  use  the  de¬ 
finition,  equation  (18),  for  negative  values  of  x ,  as  well  as  for  positive  values.  If  that  is  done, 
then  the  odd-dimensional  error  functions  tum  out  to  be  odd  functions,  while  the  even-dimensional 
functions  are  even. 


i.e., 


erf2m(x)  =  erf2n,(-x) 


(37) 


[erf2mM(x>  =  -erf2mn<-x> 
if  x  is  permitted  to  become  negative.  For  negative 
These  error  functions  are  plotted  in  figure  1  for 


x ,  the  inequalities 
positive  arguments 


(36)  will  all  be  reversed, 
and  dimensions  up  to  10. 
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Figure  1.  Generalized  Error  Rjnctionf,  erfn  (x). 


APPLICATION  TO  HYPERELLIPSOIDAL  DISTRIBUTIONS 


Retaining  the  aastunptioni  of  zero  bias  and  zero  correlations,  but  allowing  the  standard  devia¬ 
tions  to  take  on  the  unequal  vaiuea  a t,  av  .  .  . ,  an,  equation  (2)  generalizes  to 


P  (*t. 


(38) 


The  loci  along  which  p  (xj,  .  ,  .  ,  xn)  is  constant  are  given  by 


t 

4 


a  constant. 


(39) 


These  loci  thus  form  a  family  of  hyperellipsoids,  centered  at  the  origin,  with  semiaxes  equal 
to  uoj,  uo2,  ....  uan,  and  consequently  the  distribution  given  by  equation  (38)  may  be  thought 
of  as  a  hyperellipsoidal  distribution. 


Consider  now  the  probability  that  a  random  vector  r  =  (x  j,  .  .  .  ,  xn)  lies  within  the  hyper¬ 
ellipsoid 


(40) 


Prob 


X-<- 

i  ai 


"//  /  p(? 


)  dV 


hyp«r*lUpaoid 

where  dV  is  an  element  of  hypervolume  in  the  n-dimensional  space. 


(41) 


(2it)2  ox  .  .  .  on 

hyp  trail  ip  ao  id 


►  dx_ 


■// 

J  .  .  .  (2rr)l 

hyparaphara  = 

i 


dXf 


d\, 


where 


n 


(42) 


(43) 


(44) 


V 


i 


i 


t 


10 


I 


(45) 


*  * 


» 


t 


> 


introducing  hypempherical  polar  coordinates,  as  before. 
Again,  remembering  that 


Prob 


__  v  2 

y-o 


1  - 


1 


(2n)T  0 


/“  „2 

R^e’T  dR 


//  Iaa~ 


(46) 


it  follows  that 


Prob 


K 


< /< 


/ 


Rn-t  e*  T  dR 


l 


Rn_1  e“  T  dR 


P 

|VTun-1e-“2  du 

—w - 

J  u"*1  e*“!  du 


i.e., 


Prob 


<  ft1 


*  tin 


(47) 


Equation  (47)  is  the  desired  generalization  of  equation  (19)  to  the  case  in  which  the  compo¬ 
nents  of  r  have  unequal  standard  deviations. 

It  is  to  be  noted  that  the  dimensionless  quantity  /)  appearing  in  equation  (47)  reduces  to  the 

quantity  -  of  equation  (19)  for  the  special  case  in  which  “  °2  '  •  •  •  r  (,n  ° 

a 


RECAPITULATION 

The  error  function,  which  occurs  in  much  of  the  literature  of  mathematics,  physics,  and  engi¬ 
neering,  has  been  generalized  to  handle  the  normal  probability  distribution  in  n  dimensions. 

Explicit  integral  representations  for  these  functions  were  found  to  be  of  two  general  forms, 
depending  upon  whether  n  is  even  or  odd. 

Some  readily  established  recursion  formulas  and  other  relationships  were  derived  for  these 
functions. 
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